Optical response of a nematic liquid crystal cell at splay— bend 
transition: model and dynamic simulation 

Peizhi Xu/'B Vladimir Chigrinov/'0 and Alexei D. Kiselev^'H 

^Hong Kong University of Science and Technology, 
Clear Water Bay, Kowloon, Hong Kong 
^Chernigov State Technological University, 
Shevchenko Street 95, 14027 Chernigov, Ukraine 
(Dated: February 2, 2008) 

Abstract 

We study dynamical optical response of a nematic liquid crystal (NLC) cell that undergoes the 
splay-bend transition after applying the voltage across the cell. We formulate a simplified model 
that takes into account both the flexoelectric coupling and the surface rotational viscosity. The 
dynamic equations of the model were solved numerically to describe temporal evolution of the 
director profile and the transmittance. We evaluate the response time as a function of a number 
of parameters characterising dielectric and elastic anisotropics, asymmetry of the surface pretilt 
angles, anchoring energy, surface rotational viscosity and flexoelectricity. 

PACS numbers: 61.30.Gd, 78.66. Qn, 42.70.Gi 
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I. INTRODUCTION 

As it was originally shown by Berreman 
and Heffner in 1981 P|, a nematic liquid crys- 
tal (NLC) cell can be prepared to have two 
metastable states that can be switched either 
way by applying an electric field. This gen- 
eral idea underlies the mode of operation of 
bistable liquid crystal devices that have been 
attracted considerable attention over the past 
few decades. 

The approach pioneered in [l[ is based on 
using bistable twisted NLC cells that have 
two metastable twist states produced 
suit of a mismatch between the NLC pitch 
and the twist imposed by the boundary con- 
ditions at the substrates. This approach has 
been extensively studied and is found to have 
difficulties caused by fast relaxation of the 
metastable states to the intermediate stable 
configuration 3 li S, . 



An alternative approach is to use the so- 
called optically compensated bend NLC cells 

Bound- 
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also known as tt cells |7|, |8|, |9|, |1C 

ary surfaces of such cells both favour a uni- 
formly tilted alignment and the pretilt angles 
at the substrates are equal in magnitude but 
opposite in sign. For sufficiently large surface 
pretilt angles, the equilibrium orientational 
structures are non- twisted 0, S El and 
there are two director configurations that un- 
der certain conditions are degenerate in en- 
ergy: the splay (horizontal) state and the 
bend (vertical) state. 

By contrast to the bistable twisted cells, 
these bistable states are topologically distinct 
and separated by an energy barrier. So, the 
splay and bend states are both long-term sta- 
ble. Applying the voltage across the cell it 
can be switched from the splay state to the 
bend state. 

This splay-bend transition will be of our 
primary concern. We are aimed to study the 
dynamics of a NLC cell that undergoes the 
splay-bend transition induced by an external 
electric field. 



1 



The dynamical theory of NLC systems — 
the so-called nematohydrodymanics — is 
very complicated and dynamical properties of 
bistable liquid crystal cells have not received 
a fair amount of attention yet. In recent the- 
oretical studies the dynamics of vr cells |l4| . 
zenithally bistable '15^ and super-twisted jly] 
NLC devices was investigated using different 
simplified models. 

In this paper we concentrate on optical re- 
sponse of the NLC cell after switching on 
the voltage. The corresponding response 
time will be studied depending on a num- 
ber of factors such as dielectric and elas- 
tic anisotropies, asymmetry of the surface 
pretilt angles, anchoring strengths, surface 
rotational viscosity and flexoelectricity. 

The paper is organised as follows. 

In Sec. iniwe formulate our model and de- 
rive a set of dynamic equations. The numeri- 
cal results are presented in Sec. IIIII Conclud- 
ing remarks are given in Sec. IIVI 



II. MODEL 

In this section we describe our model and 
derive a set of dynamic equations. Subse- 
quently, these equations will be used to sim- 
ulate the orientational dynamics of a NLC 
layer of the thickness d that undergoes the 
splay-bend transition under the action of an 
electric field. 



constrained to lie in the x-z plane: 



n 



cos I 



where 9 is the tilt angle defined as the an- 
gle between the plane of the boundary sur- 
faces and the director; and are the unit 
vectors parallel to the x-axis and the z-axis, 
respectively. 

The vectors of easy orientation at the 
lower and the upper substrates are simi- 
larly characterised by the tilt angles 0^ and 
—6u, respectively. So, the anchoring energy 
per unit area taken in the Rapini-Papoular 
form lla is 



Wr 



anch 



sm 



Wu 



where 



e\ 



2=0, d 



(2) 

and Wl {Wu) is the 



strength of anchoring at the lower (upper) 
substrate. 

We shall also need to write the bulk part 
of the free energy per unit area 



F4n,E] = Fei[n]+Fs[n,E] 



(3) 



which is a sum of the Frank elastic en- 
ergy, Fg; [n] , and the energy of interaction be- 
tween NLC molecules and the electric field, 
Fs[n, E]. 

For the director distribution ((H), using the 
standard expression for the Frank elastic en- 
ergy [l^ gives the following result: 
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A. Free energy 

The layer is sandwiched between two par- 
allel plates, z = (lower substrate) and z = d 
(upper substrate), and we assume that both 
the electric field, E, and the z-axis are nor- 
mal to the plane of the substrates. In ad- 



dition, similar to 0, 0, 17 1, we shall re- 
strict our consideration to the case in which 
the splay-bend transition does not involve 
twisted states. 

In this case, the NLC director field, n, is 



(4) 



where dot stands for the derivative with re- 
spect to z and KeiiO) = Kn cos^ 6+Ks3 sin^ 9 
is the effective angle-dependent elastic coef- 
ficient; Kii and i^33 are the splay and the 
bend elastic constants. Similarly, the direc- 
tor field (IH) can be used to derive the expres- 
sion for the electrostatic energy FE[n, E] that 
depends on the electric field: E = -E^e^. 

Assuming that the NLC Debye screening 
length is larger than the layer thickness, the 
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NLC material can be regarded as an insula- 
tor. So, we can neglect the effects caused by 
the presence of ionic charges. 

But, the flexoelectric coupling between 
NLC and the applied field cannot be gener- 
ally disregarded. This coupling is known to 
be caused by splay and bend director distor- 
tions that give rise to an average flexoelectric 
polarisation 
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characterised by the splay and bend flexoelec- 
tric coefficients, en and 633. 

This is the well-known flexoelectric ef- 
fect, first described by Meyer in 1969 [ic^ . 
which has been extensively studied over re- 
cent years. Flexoelectricity appears to be 
a very important property of NLCs which 
must be taken into account in all experiments 
that deal with inhomogeneous director orien- 
tation. 

In our case, it is not difficult Q to obtain 
the 2;-component of the flexoelectric polari- 
sation in the following form: 

P^ = g{e)e, ^(^) = e/sin^cos^, (6) 

where Cf = 611-1-633 is the flexoelectric coeffi- 
cient. So, the final result for the electrostatic 
energy is 

FE[e,E,] = - [e,,E'j2 + P,E,]dz, (7) 
Jo 

eU9) = e^{l + usm^e), (8) 



where e 



e±6ij + (ey 



- e±)ninj is the di- 
electric tensor and u = {e\\ — e_L)/e^ is the 
dielectric anisotropy parameter. 

The Maxwell equation V x E = implies 
that the electric field E = Ez{z)ez can be ex- 
pressed in terms of the scalar potential, V: 
Ez = —V. Variation of the electrostatic en- 
ergy functional ((Tj) with respect to V gives 
the well-known electrostatic constitutive re- 
lation 



SFe 
'SE, 



,Ez + Pz = Dz 



(9) 



where is the 2;-component of the electric 
displacement field that, in contrast to E^, 
does not depend on z. 

From the relation Q the displacement 
can be expressed in terms of the voltage U = 
f^Ezdz = 1/(0) - V{d) as follows 



6iin(V ■ n) 633(11 ■ V)n, (5) where 



(10) 



6/ ln(l + u sin^ 
2ue_[_ 



(11) 



The expression on the right hand side of 
Eq. (fnH) clearly indicates the flexoelectricity- 
induced voltage shift. The effects of this shift 
in optical response of hybrid aligned liquid 
crystal cells were recently studied in [2l|. 

Since the displacement does not vary 
across the layer, it is convenient to have the 
displacement Dz as an independent field and 
use the free energy G[0, Dz\ which is related 
to the energy FJ^ E^ via the Legendre trans- 
formation I221 . 123^ 

G[e,Dz] = F[e,Ez]+EzDz, (12) 

where Ez = (Dz - Pz)/tzz- 

We can now combine Eqs. (0)-® and 
Eq. ((7j) to derive the free energy G[6, Dz] in 
the following form: 



G[d,Dz 



fbdz + fs, 



Dl 



(13) 
(14) 



tzz{ey 

fs = Unci. + Dz{ij{9o)~^{9,)), (15) 

where K{e) = K^i{6) + g'^{e)/tzz{0) is the ef- 
fective elastic coefficient renormalised by the 
flexoelectricity. 

As it can be seen from Eqs. (fT ^ -(fT ^ . the 
bulk elastic coefficient and the anchoring en- 
ergy are both renormalised by the flexoelec- 
tricity: Kei — > K and /anch fs- Static 
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properties of NLC layers submitted to an 
electric field are known to be affected by this 
renormalisation [H H 111, HQ • 



B. Dynamic equations 

Low-frequency dynamical properties of 
NLCs are generally characterised by orien- 
tational relaxation as well as by shear and 
compressional flow. A full set of dynamic 
equations governing nematohydrodynamics 
is known as the Ericksen-Leslie equations and 
describes temporal evolution of the fluid ve- 
locity and the director field. 

When the characteristic time scale of the 
velocity field is much shorter than the typical 
time of director reorientation, the flow veloc- 
ity can be adiabatically eliminated from the 
dynamics of NLC. In this approximation, the 
orientational dynamics is purely relaxational 
and can be formulated a s a time- dependent 
Ginzburg-Landau model 
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We shall apply this model to obtain the 
dynamic equation governing the orientational 
relaxation of the tilt angle in the bulk. Using 
the free energy (fT3|l gives the following result 



(16) 



where 7^ is the bulk rotational viscosity and 
prime stands for derivatives with respect to 

e. 

It should be stressed that under certain 
circumstances the backflow effect caused by 
the coupling between the fluid flow and the 
director may considerably affect dynamical 
characteristics of NLC cells. Specifically, the 
so-called "optical bounce" in twisted cells 
manifests itself as a dip in transmission of 
normally incident light after the electric field 
is turned off [Hi IHJE^ • But in cases where 
the twisted states are of minor importance 
backflow is found to induce only quantitative 
changes in the dynamics [13 . Issj ] . 



By analogy with Eq. ()16|) we can write the 
dynamic equations for the tilt angles, Oq and 
6*1, at the lower and upper substrates as fol- 
lows (il H Ei, 



7. 



09^ 
dt 



4)'A-(ft)».-f. 



0,1, (17) 



where •js is the surface rotational viscosity, 
which is defined as the ratio of the torque 
needed to change the director orientation at 
the surface for a certain angle and the corre- 



sponding relaxation velocity [37|, |3; 



Kn (N) 


6.6x10^^2 


en (C/m) 


-0.95x10^1^ 


K33/Kn 


3.0 


633 (C/m) 


-1.35x10"^^ 


d (/im) 


2.5 


e± 


6.3 


W (J/m2) 


4.0x10"^ 




12.6 


Ol (deg) 


46.0 


uiy) 


20.0 


Ou (deg) 


44.0 


no 


1.5 


7b (N s/m2) 


0.1 




1.6 


75 (N s/m) 


3.0x10-^ 


A (pm) 


0.55 



TABLE I: Parameters of the model employed in 
the calculations. 



III. SIMULATION RESULTS 

In this section we present our numerical 
results obtained by solving the dynamic equa- 
tions (fT^ and (fT7|) numerically. Dependen- 
cies of the tilt angle on z at specified points 
in time were computed using the finite differ- 
ence time domain method. The parameters 
used in our calculations are listed in Table IB 

In order to study the dynamics of opti- 
cal response of the layer, the data represent- 
ing temporal evolution of the director pro- 
file, which is the tilt angle function of 
9{z,t), were used as an input for computing 
the transmittance of light through the layer 
placed between two crossed polarisers. 

The expression for the transmittance 
can be derived by using the Jones matrix 
method 13911. When the director of LC cell 
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FIG. 1: The director configuration through the cell at different points in time after applying the 
voltage. The anchoring strengths at the substrates are assumed to be equal, Wl = Wu = W, and 
the list of the parameters are given in Table H 



is at 45 degrees to the input polariser, the 
transmittance, T, is given by (41i, |41|J 



T = sin2(A0/2), 
1 sin^ 9 cos^ 9 



A(f) 



n 



eff 



nt 



nt 



(18) 
(19) 

(20) 



where A0 is the phase difference between 
the ordinary and extraordinary ray; A is the 
wavelength of the incident light and no {rie) is 
the ordinary (extraordinary) refractive index. 

Thus, we describe the dynamics of opti- 
cal response by computing temporal change 
in the transmittance (jl8p . An important pa- 
rameter characterising the rate of change of 
the transmittance is the response time which 
is the time it takes for the transmittance to 
increase from 10% to 90%. 



We begin with the case in which the flexo- 
electric effect is neglected and Cf = 0. Fig. [T] 
shows how the director profile evolves in time 
after applying the voltage across the NLC 
cell. Asymmetric pretilt angles and symmet- 
ric surface anchoring energy are used in this 
calculation. As is illustrated in Fig. ^ the 
initial director configuration corresponds to 
the splay state which gradually transforms 
into the bend state under the action of the 
electric field. 

Now we pass on to discussing the ef- 
fects related to the dielectric and elastic 
anisotropics. The results for various values 
of the dielectric anisotropy parameter, u = 
(^11 ~ e±)/e±! and the elastic ratio, K^^/Ku, 
are shown in Figs. El and El respectively. 

Fig. |2(b)| indicates that the response time 
is a non-monotonic function of the dielec- 
tric anisotropy parameter and goes through a 
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FIG. 2: (a) Transmittance as a function of time at various values of the dielectric anisotropy 
parameter, u = (en — e_L)/e_L. (b) Response time as a function of u. 



minimum in the vicinity of m = 0.8. By con- the substrates (fTTj) . 
trast, as is shown in Fig. |3(b)[ the response 
time monotonically dechnes as the ratio of 
i^33 and Kii increases. So, large values of the 
elastic ratio facilitate the splay-bend transi- 
tion. 



The surface pretilt angles, 9l and 9u, are 
known to play an important part in the splay- 
bend transition [s^ . These are among the pa- 
rameters that affect the dynamics of optical 
response through the boundary conditions at 



The first parameter we consider is the dif- 
ference between the pretilt angles: A^^ = 
9l — du- Fig- |4(a)| shows the curves for the 
transmittance varying in time at various val- 
ues of the pretilt angle difference. As is ev- 
ident from Fig. |4(a)[ the curves are getting 
steeper as AO^ increases and the response 
time, shown in Fig. |4(b) is a decreasing func- 
tion of AOs- 

The anchoring energy dependence of the 
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FIG. 3: (a) Transmittance as a function of time at various values of the elastic ratio A'sa/i^n, 
(= (dh — (^u))- (b) Response time as a function of the elastic ratio K33/K11. 



response time is plotted in Fig.^lfor the sym- 
metric case with = Wu = W. The 
curve is depicted in logarithmic scale and 
clearly indicates the transition between two 
regimes of anchoring: the weak anchoring 
regime and the strong anchoring regime. In 
the regime of weak anchoring, the extrapola- 
tion length is larger than the cell thickness, 
d, and the response time is small. As is seen 
from Fig. ^ the response time increases with 
the anchoring energy and saturates on reach- 
ing the strong anchoring regime where the 



extrapolation length is much smaller than d. 

Influence of asymmetry in the anchoring 
energy strengths on the response time is illus- 
trated in Fig. iniwhere the anchoring strength 
at the upper substrate is kept constant at 
the value listed in Table ^ Wu = W. It is 
shown that the response time varies slowly 
and reaches its maximum at Wl/Wu ^4.0. 

The surface rotational viscosity, 7^, can be 
conveniently characterised by the ratio of 7^ 
and the bulk viscosity which has the dimen- 
sion of length. There are, however, only few 
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FIG. 4: (a) Transmittance as a function of time at various values of A^^, (= 9l — 9u). (b) Response 
time as a function of A^^. 



measurements of this length that, according 
to [s^, 13]) can be of the order tens and 
hundreds nanometers. Our numerical results 
on the surface viscosity dependence of the re- crystals j4 
sponse time are presented in Fig. [7| It is seen 
that variations of the surface viscosity over a 
wide range of values have almost no effect on 
the response time. 



appears to be eliminated from the consider- 
ation. There are some measurements of the 
flexoelectric coefficient in a variety of liquid 



SHHHIii • It was found 
that the value of |e/| typically falls in the 
range between 5 x 10~^^ C/m and 9 x 10~^^ 



So far we have limited our discussion to 
the case in which the flexoelectric coefficient 
6/ vanishes and thus the flexoelectric effect 



C/m. But reliable and accurate experimental 
estimates of ej are still missing. For example, 
the reported values of for MBBA turned 
out to differ in both magnitude and sign de- 
pending on theoretical approach used for pro- 
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FIG. 5: Response time as a function of the anchoring strength, W = Wl = Wjj. 



cessing experimental data j21, 2^ isl. l45l|. 

Numerical results related to the effect of 
fiexoelectricity on the dynamics of NLC cell 
are presented in Figs. I5HTU1 The curves 
shown in Fig. |H1 indicate that the dielectric 
anisotropy dependence of the response time 
turns out to be strongly affected by the fiex- 
oelectric effect. In the presence of fiexoelec- 
tricity the curve has a pronounced maximum 
peaked at u ~ 0.5 which follows a minimum 
reached at m ~ 0.35. 

By contrast to the dielectric anisotropy 
dependence, the dependencies of the response 
time on the pretilt angle difference, A6's, de- 
picted in Fig. ini do not differ significantly. 
For AOs larger than 10 deg, referring to 
Fig. ini the curve with non-zero fiexoelectric 
coefficient is approximately shifted upward 
by 5 ms with respect to the curve computed 
at vanishing e/. 

Finally, we comment on the dependen- 



cies displayed in Fig. UHl The curves plot- 
ted in Fig. |10(a)| represent temporal evolu- 
tion of the transmittance at different values 
of the fiexoelectric coefficient. The response 
time in relation to the fiexoelectric coeffi- 
cient obtained from these curves is shown in 
Fig. |10(b)} It can be seen that the response 
time steeply declines after reaching a maxi- 
mum at 6/ ^ —1.15 X 10"^^C/m. 



IV. CONCLUSION 

In this paper we used a simplified ap- 
proach to study the dynamics of optical re- 
sponse at the splay-bend transition that oc- 
curs after applying the voltage across the 
NLC cell. It is assumed that the coupling be- 
tween the director and the flow velocity can 
be eliminated from consideration. 

Similar approach was recently applied 
to formulate the model of switching in 
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FIG. 6: Response time as a function of the anchoring strength ratio, Wl/Wjj, at Wu 
4.0 X lO-^J/m^. 
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a zenithally bistable device [15[. In our 
case, however, not only the boundary condi- 
tions (fT7|) are different, but also inhomogene- 
ity of the electric field is taken into account 
using the constitutive relation 0. 

The simulation results for the transmit- 
tance were obtained by solving the dynamic 
equations of the model numerically. The re- 
sponse time characterising the rate of change 
of the transmittance was evaluated to study 
how the parameters of the cell influence the 
dynamics of optical response. 

Dependencies of the response time on 
the dielectric anisotropy parameter and on 
the fiexoelectric coefficient are found to be 



strongly non-monotonic. It was shown that 
the response time declines as the elastic ra- 
tio K^2,lKii or the pretilt angle difference 
AOs increases. From the other hand, the re- 
sponse time appears to be relatively insensi- 
tive to anchoring strength asymmetry and to 
changes in the surface viscosity. 
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